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We propose the formulation of lattice QCD wherein all elements of the theory (gauge 
action, fermionic action, 9- term, and all operators) are constructed from a single object, 
namely the lattice Dirac operator D with exact chiral symmetry. Several regulariza- 
tions of this type are suggested via constructing scalar densities (gauge actions) that 
are explicit functions of D. The simplest of these is based on the proposition that clas- 
sical limit of gauge density associated with trace of D is (up to an additive constant) 
proportional to tr F 2 , while the corresponding operator is local. The possibilities of ex- 
plicit interrelations between gauge and fermionic aspects of the theory are emphasized 
together with the utility of such formulations for exploring the QCD vacuum structure. 

o 

1 Introduction 

Ph. 

This is the second in the series of papers whose purpose is to introduce the set of approaches 
aimed at creating a systematic framework for studying QCD vacuum structure in the path 
integral formalism [1]. The basic premise of such framework is the use of lattice QCD both to 
define the relevant representation of the vacuum (the associated ensemble of Euclidean QCD), 
as well as an exclusive source of necessary dynamical information obtained via numerical 
simulation (the Bottom-Up approach). One of the important points emphasized in Ref. [1] 
is that in the search for the fundamental structure (the structure relevant for all aspects 
of QCD physics), it is beneficial to utilize the freedoms we have in the lattice definition 
of QCD. In particular, one can take advantage of the fact that the degree of space-time 
order in typical configurations can vary significantly over the set of valid lattice theories at 
a given cutoff. 1 A possible tool for reaching theories with high level of space-time order is 
the use of particular transformations in the set of gauge configurations, chiral orderings [1], 
which replace given link U n ^ with the effective SU(3) phase U n ,n associated with hopping of 
properly defined fermion from site n + fi to site n. More precisely, it was suggested that such 
configuration-based deformation of the gauge action tends to preserve the physical content 
of the lattice theory [principle of chiral ordering) while increasing the space-time order. 



1 A "degree of space-time order" in a given configuration can be defined in principle via Kolmogorov 
entropy (algorithmic complexity) of binary strings representing its coarse-grained descriptions [1]. 
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The strategy of searching in the set of actions via chiral ordering transformations takes 
advantage of the fact that, for purposes of the space-time structure, it is not necessary 
to know the form of the underlying action explicitly. Indeed, all we need is the access to 
its ensemble ("typical configurations"). However, it would certainly be beneficial if it was 
possible to work with actions that generate a high degree of space-time order and, at the 
same time, can be explicitly written down. In this paper we suggest the class of actions 
that might satisfy this requirement. To do so, we will adopt a characteristic feature of the 
chiral ordering approach, namely that the transformed gauge actions become functions of 
the lattice Dirac kernel D on which the transformation is based. In other words, we propose 
gauge actions that are explicit functions of D. While this does not in itself guarantee that 
the typical configurations in these theories will be significantly more ordered than those of 
standard actions (e.g. Wilson gauge action), we think this is a very reasonable expectation. 
Moreover, the theories so obtained will facilitate various interrelations between gauge and 
fermionic aspects of the theory and are interesting in their own right. We emphasize (as we 
did in Ref. [1]) that our approach acquires its full potential only in the framework of QCD 
with dynamical fermions. Definition of full theory based on a single object, namely a lattice 
Dirac kernel D describing chiral fermion, is perhaps a cleanest manifestation of this fact. In 
what follows, we will refer to any lattice regularization based on this idea as coherent lattice 
QCD. Also, "lattice QCD" will frequently be abbreviated as LQCD. 

Upon accepting the possibility of coherent LQCD, one realizes rather quickly that there is 
a large freedom of choice available here. We explore this freedom to some extent by searching 
for a formulation that casts gauge and fermionic contributions into an overall dynamics of 
the theory in the most form-symmetric manner. In the resulting regularization, symmetric 
logarithmic LQCD, the gluonic contribution can be viewed exactly as that of an infinitely 
heavy quark. Conversely, the full effective action can be expressed as a sum of Nf + 1 
gluon-like terms with only one being local (representing the usual gluon), and the rest of 
them non-local. The meaning of non-local F 2 is that for smooth configurations the operator 
effectively averages F 2 over the physical distance associated with quark mass in question. 
The feasibility of symmetric logarithmic LQCD rests on the resolution of certain locality 
issues and the ways of approaching them are discussed in detail. 

One can alternatively arrive at coherent LQCD via formal/classical considerations in the 
continuum. Indeed, it can be argued that the full action for a classical configuration (smooth 
almost everywhere) can be written entirely in terms of Dirac operator D. The particular 
way in which the corresponding formal equation comes about suggests a prescription for 
lattice regularization of QCD once lattice Dirac operator is specified. We refer to this class 
of lattice theories as classically coherent LQCD due to the important role of classical fields 
in obtaining it. 

In order to make the lattice theory fully coherent, it is desirable that the operators used 
for measuring physical observables are also explicit functions of D. This has, in fact, been 
already achieved in a generic way via the use of a chiral ordering transformation [1]. However, 
it would be useful if at least certain important local composite fields could be expressed in 
yet more explicit way. We discuss in some detail the possibilities offered by considering 
various Clifford components of D. 

One of the interesting aspects of coherent LQCD is that it admits a natural definition of 
effective LQCD at fermionic response scale along the lines discussed in [1]. Let us recall 
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that the purpose of effective LQCD is to study the effective structure of QCD vacuum, where 
the influence of fluctuations above a certain scale is suppressed. Effective structure defines 
an "unfolding" of all-encompassing fundamental structure into the one that is relevant for 
processes that "excite" mainly fluctuations down to some desired length scale. Effective 
LQCD at fixed fermionic response scale is the first step in such a definition. The feature 
of coherent LQCD is that, contrary to the generic case, it allows for definition of effective 
LQCD in an explicit manner, i.e. via the action that can be written down and, at least in 
principle, directly simulated. This is discussed in the last part of the paper. 

2 Coherent Lattice QCD 

In this section we discuss the simplest version of coherent LQCD, where the gauge action 
is constructed in direct analogy with representation of topological term using lattice Dirac 
operator with exact chiral symmetry. As a basis for such construction we propose the 
following conjecture. 2 

Conjecture C3. Let A^x) be arbitrary smooth su(3) gauge potentials on R 4 . IfU(a) = 
{ U n ^{a) } is the transcription of this field to the hypercubic lattice with classical lattice 
spacing a, and I = { U n ^ — > I c } is the free field configuration then 

tr (D ofi (U(a)) - A),o(I)) = ~c s a 4 trF^(0)F^(0) + 0(a 6 ) (1) 

for generic D e S F . Here c s is a non-zero constant independent of A^{x) at fixed D, and 
F^ v {x) = d^A u (x) — d u A^(x) + [A^(x), A v (x) } is the field-strength tensor. 

We wish to make the following remarks for completeness. 3 

(i) By potentials being "smooth" we mean differentiable arbitrarily many times. Relation 
(f) is expected to be valid also for configurations with singular gauge potentials A^(x) if 
they are smooth in some finite neighborhood of x = 0. 

(ii) Transcription to the hypercubic lattice with classical lattice spacing a is defined via 

U n A a ) = Vexp^aJ^ ds A^(an + (f - s)a/t)) (2) 

where V is the path ordering symbol and ft a unit vector in direction /i. 

(Hi) The set S F , defined in Ref. [1], contains single-flavor lattice Dirac operators D (or 
equivalently actions ipDip) satisfying the usual requirements including the exact lattice chiral 
symmetry. 

2 We will keep the notation consistent with the first paper of the series (Rcf. [1]) to the largest extent 
possible. The numbering of conjectures and definitions is also continued. Conventions for continuum gauge 
theory are summarized in Appendix A for convenience. 

3 It is desirable to specify and fix our convention for denoting traces. The basic rule is that tr denotes 
a local trace, i.e. trace of a matrix at given x or n, while Tr denotes a global trace, i.e. trace involving 
space-time coordinates. In the above form, the traces are implicitly assumed to be taken in the linear space 
"natural" to the object in question, e.g. in Eq. (1) the trace on left-hand side is over spin-color, while the 
trace on the right-hand side is over color only. If it is still necessary to distinguish spin and color traces, we 
use the corresponding superscripts, e.g. tr s , tr c . 
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(iv) Due to the translation invariance of D it is sufficient to consider the statement of 
Conjecture C3 only at x = n = 0. Analogous statement is true for arbitrary fixed x. 

(v) The underlying reason for expected validity of Conjecture C3 is that trD nn is scalar 
(under hypercubic group), local, gauge invariant function of the gauge field. Thus, the 
continuum operators appearing in its asymptotic expansion in classical lattice spacing will 
be the gauge invariant operators (of the appropriate dimension) which are also scalar under 
hypercubic group. Up to dimension four, the possibilities for such operators in the continuum 
only include a constant and tr F pv F pu . This leads to the proposed classical limit. 

(vi) The relation analogous to (1) for the pseudoscalar case [2] (see also Ref. [3]), namely 

tr 75 (L> ,o(C/(a)) - D 0fi (lj) = -c p a A trF^(0)F^(0) + 0(a e ) (3) 

is a basis for definition of topological charge density via chirally symmetric Dirac operator. 
Here F pv = \e pupa F p(T (with ei 2 34 = 1) is the dual of the field-strength tensor, and the 
constant term on the left-hand side vanishes. 

(vii) For the family of standard overlap Dirac operators [4], based on Wilson-Dirac operator 
with mass — p, the existence of constant c p (p) has been verified in explicit calculations [5] 
and is given by 4 

^ = A 0<„<2 (4) 

The validity of Conjecture C3 for the family of overlap Dirac operators will be studied 
explicitly in Ref. [6]. The corresponding constants c s (p) will be determined there using both 
analytical and numerical methods. 

(viii) In practice, we frequently consider QCD in finite physical volume and the statement 
analogous to Conjecture C3 can be verified most directly in such setting. The relevant 
conclusion can be formulated as follows. 

Conjecture C3a. Let A p (x) be arbitrary smooth su(3) gauge potentials on symmetric torus 
of size L p , and let = L p /L be the classical lattice spacing on L A hypercubic lattice. If 
U(a L ) G U L is the transcription of A^(x) to this lattice then the following non-zero finite 
limit exists 

tr(D 0fl (U(a L )) - DqM) 

ai tr iV(0)iV(0) ~ [b) 
for generic D G S F . Here c s is a constant independent of A^x) at fixed D. 

Due to the locality of D, the constants c s in Conjectures C3 and C3a are expected to be 
identical. 

We can now define the action for simplest version of coherent LQCD with Nf flavors of 
quarks and the CP-violating ^-term as 

N f 

S mmf} = Tr(P-ie l5 )(D(U)-D(I)) + Y,^{ D (U)+m f )^ (6) 

/=i 



Note that the "standard" overlap operator is given by £)W = p[l + A(A^A)-i] where A = D w - p and 
Dw is the massless Wilson-Dirac operator. 
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where {rrif} is the set of real non-negative quark masses, and 5 



m = T^ = ff W=iSP? ("e (-»,»]) (7) 



Indeed, the locality of S 1 defined above follows from locality of .D, the required symmetries 
follow from its transformation properties, and the classical correspondence is equivalent to 
Eqs. (1) and (3). Note that the free-field term in the gauge part of the action only contributes 
a field-independent constant and can be discarded. 

The dynamics of coherent QCD (6) is completely encoded in the chirally symmetric 
lattice Dirac operator. This feature appears particularly clearly after fermionic variables are 
integrated out, which leads to the distribution density of the gauge fields given by 

PpAm oc e ^fHD+m) + (-p+ih5)D] = det [(D + m) N f e^ +iS ^ D } (8) 

where D = D{U), and we consider lattice QCD with Nf mass-degenerate flavors of mass m. 
This can also be written as 



P Mm oc det[e^ D ^(-±) Nf e-^ 



det 



d/3 

-j3(D+m) I " \ N f 



{-^) ^ {D+m) ] ( 9 ) 



which emphasizes the close explicit relation between fermionic and gauge contributions to 
the distribution density of full QCD. 

Let us finally comment on three points. 

(a) Note that the gauge action density proportional to tr(D(U) — D(l)) n ^ n is expected to 
be a non-ultralocal function of U. Thus, in coherent QCD we are necessarily giving up the 
possibility of strict reflection positivity at non-zero cutoff. However, this feature was aban- 
doned already by requiring exact chiral symmetry of the fermionic action since such actions 
are necessarily non-ultralocal in fermionic variables [8]. Reflection positivity is expected to 
be recovered without problems in the continuum limit. 

(($) Due to a complicated nature of Ginsparg-Wilson operators, it is expected that there 
exist configurations U for which D(U) is not uniquely defined. This can happen, for example, 
in fine-tuned backgrounds, where the number of zeromodes (topology) changes. For such 
configurations the coherent gauge action is also not expected to be well-defined. We thus 
have to take these configurations out of the path integral (assign a zero distribution density 
to them) and operate under the assumption that we still obtain a universal theory (QCD) 
in the continuum limit. Again, however, we have made that choice already by using chiral 
fermions to begin with. 

(7) A noteworthy aspect of coherent LQCD is that the formulation puts the fermionic and 
gauge aspects of a numerical simulation with dynamical fermions on the same footing. The 
fact that fermionic and gauge parts of the action respond more coherently to proposed 



5 We should perhaps emphasize here that everything in this equation, including the bare parameters, is 
in lattice units. See also the notational conventions set in Ref. [1]. 



5 



changes of the configuration in the Markov chain could be beneficial for algorithms that may 
be used in their simulation. Moreover, while adopting Ginsparg-Wilson fermions certainly 
makes the problem of dynamical simulations more involved, putting the gauge part of the ac- 
tion on par with fermions doesn't appear to generate qualitatively new complexities. Indeed, 
the simplest coherent LQCD discussed in this section requires simulating det(D + m) N f e _/3D 
instead of det(D + m) Nf , and this can be incorporated naturally in the existing algorithms. 
Detailed discussion of related issues will be given in Ref. [7]. 



3 More General Coherent LQCD 

While the form of coherent LQCD introduced in Sec. 2 is perhaps the simplest one, there 
is in principle a much larger set of possibilities here. Indeed, by coherent LQCD we mean 
any formulation where the gauge and fermionic parts of the action are tied together in an 
explicit manner. Since the Dirac operator D defining fermionic action is a much more general 
(and yet highly constrained) object than the structure associated with the gauge action, it 
is natural to view D as primary and to model the gauge part accordingly. 

Following this route, we can obtain a large class of coherent LQCD formulations if we 
treat general functions f(D) in a similar manner as we treated D in Sec. 2. In particular, 
if f(D) is a local operator and if tr/(D) n „ is a scalar lattice field, then we generically 
expect the validity of a statement analogous to Conjecture C3 with equation (1) for smooth 
backgrounds replaced by 



tr 



f(D(U(a))) Qo - /(£>(!)) ] = -cV tr ^(0)^(0) + 0(a 6 ) (10) 



with the appropriate constant c s . Proceeding in the same way also for the pseudoscalar case 
(i.e. considering operator 75/(1}) ) we obtain the definition of coherent LQCD in this case 
by replacing D in the gauge part of Eqs. (6,8) with f(D), while Eq. (7) remains unchanged. 

The freedom in choosing f(D) may allow us to formulate lattice QCD in a manner that 
casts the gauge and fermionic parts of the full action into an even more mutually symmetric 
form. To do that, let us consider the function f(D) = ln(D + with fixed r\ > 0, and where 
the principal branch of the complex logarithm is used to define the operator. The analog of 
Eq. (6) is then (up to an inessential constant) 

S mmf} = Tr(P-ie l5 )\n(D(U)+m ) + £ &(d(U) + m f )^ (11) 

/=i 

where we set rj = m > to reflect its mass-like nature but still put it in contrast to quark 
masses rrif. Eq. (8) for distribution density with degenerate quark masses is replaced by 



Psb CX e r ^i N fHD+m) + (-0+i§ 1 5)HD+m o )} = ^ 



(D + mff (D + mo)-^" 75 (12) 



where D = D(U) and the formal relation A B = e^ nA ^ B for square matrices A and B was 
used. We will refer to the above regularization as logarithmic LQCD for obvious reasons. 

The validity of logarithmic LQCD relies on the following two conjectures (classical limit 
and locality) that we propose for the numerical and analytical investigation. 
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Conjecture C4. Let A^(x) be arbitrary smooth su(3) gauge potentials on M 4 . IfU(a) = 
{ U n ^{a) } is the transcription of this field to the hypercubic lattice with classical lattice 
spacing a, and 77 > 0, then equation (10) is valid for f(D) = \n(D + 77) with generic D e S F 
such that f(D) is well-defined. The non-zero constant c s = c s (r]) is independent of A^(x) 
at fixed D . 

Conjecture C5. Let 77 > and D e 5 F such that f(D) = \ia(D + 77) is well-defined. Then 
s n (U) = tr f(D(U)) n>n is a local composite of U . 

By saying that f(D) is "well-defined" in the above conjectures we mean that f(D{U)) 
is a uniquely defined linear operator for all backgrounds U for which D{U) is uniquely 
defined. This is true e.g. for the family of overlap Dirac kernels D^ p \ It should be also 
mentioned that by local composite field s n (U) we implicitly understand the strong version 
of locality condition which can be formulated as follows. Representing the configuration U 
in the canonical form used in Ref. [1], where U ntfJi <-> u{n, ji) = {u a (n,fi), a = 1,2,..., 8} 
with u a (n, ji) being independent real parameters, locality of s n (U) requires the existence of 
a(rj) > and A(rj) > (independent of U) such that 



max 

fi,a 



ds n {U) 



du a (m,n) 



< A(ri)e~ aM \ n ~ m \ Wn,m (13) 



and for all configurations U for which D(U) is uniquely defined. 6 It needs to be emphasized 
that locality in this form has not been proved rigorously even for simpler cases such as tr D$ n 
with being an overlap Dirac operator. However, one can at least support it by exploring 
the behavior numerically in some average sense or for a specific class of configurations as was 
done e.g. in Refs. [9, 10] for the case of fermionic locality. Note also that the coefficients a, 
A are implicitly assumed to be chosen in the "optimal way", which is in principle obtained 
by fixing the maximal possible a first, and then selecting the minimal A given that choice. 

The heuristic argument for the validity of Conjecture C5 is as follows. Consider the 
asymptotically large value of rj. In that case we have that \n(D + rj) = In 77 + ln(l + D /rj) 
In 77 + D/rj. Thus for r] — > 00 we expect exponential localization of tr[ln(D + r))] n>n with 
A{r]) — > and 01(77) ~~ a D, where ap is the localization range of tr D nn . For finite 77 > 0, 
the operator ln(l + D/rj) = Ylx^x^-i^ + ^/v)^ (with Dip\ = \ip\) is expected to be 
both well-defined and local since nothing non-analytic happens when performing a spectral 
transformation A — > ln(l + A/77) if the spectrum is bounded and if there are no negative real 
eigenvalues (< —77). Both of these conditions are generically satisfied by elements of S F . 
For example, if we consider the overlap Dirac operator with p — 1, one would naively expect 
that something non-analytic could happen for 77 < 2 because the expansion of ln(l +x) does 
not converge for \x\ > 1. However, this is not the case (we do not rely on the expansion 
at all) since there are no singularities of the complex logarithm in the relevant part of the 
complex plane. The situation in this regard is similar to (D + r/)^ 1 = r/^ 1 (1 + D / r/)^ 1 which is 



6 Note that this definition implicitly assumes that we can exclude, without any harm, configurations for 
which D(U) is not uniquely defined. In addition, one might have to exclude a subset of configurations for 
which the condition (13) is not satisfied if it can be shown that these configurations will not contribute in 
the path integral. This makes both the rigor and the possible numerical verification of strong locality rather 
difficult to deal with if non-ultralocal operators are used. 
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expected to be exponentially localized for arbitrary r\ > (exponential decay of the massive 
propagator) and an analogous argument applies. 

We now wish to make the following remarks: 

(%) Note that for fixed D £ S F the mass-like parameter m merely labels different formu- 
lations of coherent LQCD with f(D) = ln(D + mo), and is fixed as the continuum limit 
is taken. The freedom in choosing m can be used to adjust the desired degree of lattice 
locality in gauge interactions. 

(ii) We have implicitly assumed so far in our considerations that both parts of gauge inter- 
actions (scalar action and the #-term) are governed by the same function f(D) . This is 
natural from the point of view of "coherence" between different parts of the theory. How- 
ever, it might prove useful in some circumstances to use different functions f s (D) and f p (D) 
for the scalar and pseudoscalar parts respectively. For example, in the case of logarithmic 
LQCD we would have 



Pp,e,m oc det 



(D + m) N f (D + mf (D + m p f^ 1 (14) 



instead of (12). We emphasize that the distributions in question are actually invariant under 
the change of f p (D) (the global topological charge will not change), but the local behavior 
of the action will be different, which might be useful for certain theoretical considerations. 

(Hi) While the relation between the scalar constant c s (D) corresponding to D, and c s (f(D)) 
is in general complicated, the pseudoscalar case can be very simple. For example, if D is an 
overlap operator with p — 1, then the pseudoscalar constant of tr 75/(1?) is given by 

c P U\D)) = f{2) ~ f{0) c p {D) (15) 

where we have assumed that [D, f(D)] = 0. The generalization to arbitrary D £ S F in 
terms of its real modes is straightforward. We emphasize that the strict topological nature 
of ti^Dn^n is inherited in the tr 7s/(I?)n,n- 

(iv) Another choice of function f(D) that could be used to obtain coherent LQCD with 
tunable range of lattice locality is f(D) — (D + r])^ 1 . This possibility, while appearing 
impractical, has an interesting theoretical appeal since it involves an operator (quark prop- 
agator) which enters expressions for hadronic correlation functions. 

(v) Conjecture C4 offers a straightforward insight into the nature of fermionic determinants 
for classical backgrounds. This is discussed in Appendix B. 

(vi) It is amusing to note that, as can be seen from Eq. (12), the role of continuum limit 
driving parameter (3 in logarithmic LQCD is analogous to that of the number of quark 
flavors. Indeed, (3 can be viewed as counting the (continuous) number of additional fermions 
(if c 5 (m ) < 0) or pseudofermions (if c s (m ) > 0). The process of taking the continuum 
limit |/3 1 — > oo involves engaging more and more of such particles. In the continuum limit 
there is infinitely many of them and they become infinitely heavy since mo is kept fixed and 
thus mo/a — > oo as a — > 0. 
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(vii) Finally, we note that the straightforward generalization of Conjecture C3 is expected 
to be valid for arbitrary polynomials of D. For future reference (see Sec. 5), let us make the 
needed statement explicit. 

Conjecture C3g. Let A^x) be arbitrary smooth su(3) gauge potentials on M 4 . IfU(a) = 
{ Un,[i{a) } is the transcription of this field to the hypercubic lattice with classical lattice 
spacing a, then equation (10) is valid for f(D) = D n , where n is a positive integer and D is 
a generic element of S F . The non-zero constant c s is independent of A^(x) at fixed D. 

4 Symmetric Logarithmic LQCD 

Following the ideas outlined in the previous section, one can attempt to define coherent 
LQCD in the form which is even more symmetric with respect to fermions and bosons. This 
can be done by exploiting the expected behavior of functions c (if) and a(rj), A{rj) associated 
with s n = trln(D + rf) n , n in the i] —> limit. It turns out that the resulting formulation 
that we present below is not expected to be local by the usual criteria of strong locality. 
However, we will argue that if the locality condition is weakened to the form which still 
appears physically acceptable (weak locality), our symmetric formulation could be considered 
local, and thus expected to define the correct continuum theory at least for some minimal 
number of quark flavors. Given the elegant appeal of this regularization and the insight it 
offers into an interplay between bosonic and fermionic degrees of freedom, we analyze these 
issues in some detail below. 

Before we begin, let us recall that we have already concluded about the functions c s (rj), 
a(rj), A{rj) of s n that 



with c s D , «d, Ad denoting the corresponding constants for D G S F . Also, c s (rj), a(r)) and 
A(rf) are generically expected to be non-singular for rj > 0. As for the behavior near the 
lattice massless limit 77 — ^ 0, one can make some very plausible guesses. We begin with 
discussion of c s (rf) in this regard, which will bring us directly to the symmetric formulation. 
Since the norm of \n(D + if) will diverge logarithmically in the 77 — ^ limit for smooth 
configurations, it is expected that s n , and thus c s (i]), will diverge at least logarithmically as 
well. We thus propose the following conjecture to be verified via explicit calculations. 

Conjecture C6. Let D e S F such that f(D) = \n(D + rf) is well-defined for arbitrary 
r] > 0. If c s (i]) is the associated classical coupling of s n (U) = tr f(D(U)) n , n to tr F^F^ 
defined by Conjecture C4 7 then there exists i] such that c s (rj) is monotonic for i] < i] , and 
lim^o |c 5 (r])| = 00. Moreover, there exists a non-zero (possibly infinite) limit 



lim c s (ti) = lim — = 



lim aivj) = old 




(16) 




hm^M 
f7->o ln(r?) 




K S (0) ^ 



(17) 
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Note that the possibility k s (0) = ±00 is included here. With regard to Conjecture C6 it is 
worth mentioning that the anticipated divergence is manifestly present for the pseudoscalar 
constant c p (rf). For the family of overlap Dirac operators one can see easily that 



cP (") = i^ 111(1 + 7» kP(0) = -t^ (18) 

4.1 Definition of the Lattice Action 

The divergence in coupling of trln(D + r\) n ^ n to F 2 implied by Conjecture C6 suggests that 
we can trade the continuum-limit driving coupling (3 of formulation (11,12) in favor of 
"variable" gauge mass parameter m Q = m (g). 7 To see this explicitly, consider the action 
after integrating out fermions (effective action) of formulation (11,12) with 9 = 0. 

-S% = Tr[N f ln(D + m)-(3HD + m )] p = (19) 

The correct classical limit at arbitrary g will be reproduced when we replace it with 

S e J f mo =Tr[N f \n(D + m) - sgn(c 5 (m )) ln(£> + m ) ] (20) 
where sgn(a;) denotes the sign function, and m = mo(g) is implicitly related to g via 

\ cS ( m o)\ = ^ (21) 

Note that the Conjecture C6 implies the existence of g > such that g G (0, g ] is in 
one-to-one correspondence with m G (0, 770] . 

Let us now write this lattice regularization explicitly in its general form. One way to 
do that is to emphasize the fact that, instead of g, 9, {m^}, the parameters of QCD in our 
formulation come all naturally in the mass-like manner, namely mf, tUq, {rrij} . Indeed, the 
action of symmetric logarithmic LQCD can be defined by (up to an inessential constant) 

Sms >m p ,{m f} = sgn(c 5 K)) Tr \n(D(U) + m$) 

- i sgn(c p (|m^|)) sgn(m^) Tr 75 \n(p(U) + |m^|) 

Nf 

+ J^$ f (D(U)+m f )^ f (22) 
/=i 

where g is related to uIq via Eq. (21), while 9 is related to through 

9 = sgn(m^) 16tt 2 \c p (\m£\)\ 9 e (-tt, tt] (23) 

Note that in the above definition we have allowed the pseudoscalar mass m p to be negative 
so that the lattice action density preserves exactly the transformation property of the pseu- 
doscalar part under 9^—9. We emphasize that lin^p^^ c p (m p ) = 0, and 9 = case is 



7 We thank Andrei Alexandra for sparking this line of thought in our conversations. 
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thus obtained in this limit (as well as in m p — > — oo limit). The above definition assumes 
that |c p (7/)| is a decreasing function of 77 in the range of positive values [m^' c , 00) such that 
9(m P,c ) = 7r = 16n 2 \c p (m P,c )\. This is true for the family of overlap Dirac operators D^ p \ 
in which case c p {m p ) is given by Eq. (18) and we have 

The mass-like parameters in the above lattice theory are thus chosen to vary in the ranges 

rrif £ (0,oo) mo G(0,m ' c ] m P E[m P,c ,oo)U(—oo,—m P,c ) (25) 

where mo' c is the maximal /xq satisfying the statement of Conjecture C6. The continuum 
limit with Nf < 16 is taken via = m^{a) — > while decreasing rrif = m,f(a) oc arh^ 
towards zero so that some set of renormalized masses (in physical units) is held fixed, 
and with the pseudoscalar mass m p (specifying 6) kept constant in the process. 8 

The second possibility of a general definition that is worth mentioning is obtained by 
setting = rn p = m with the scalar part of the gauge action being identical to that of (22) 
(with continuum limit driven by changing mo according to (21)), but with the pseudoscalar 
part being constructed as in logarithmic LQCD (11,12). The virtue of such definition is that 
both scalar and pseudoscalar action densities then involve gauge operators with the same 
locality properties driven by m . This is relevant if one is interested in studying the space- 
time behavior of the action. We emphasize again that the actual distribution density P(U) 
(and thus the typical configurations) is invariant under various choices of #-term constructed 
from fixed D G S* . Indeed, all such choices yield by construction the same global topological 
charge Q(U) as the one defined by D(U). The different formulations here will differ only by 
the space-time behavior of the pseudoscalar (topological) action density. 



4.2 The Locality Issue 

Let us now turn to the expected behavior of 0(77) for 77 — > 0. Here one expects that 
lim^o 01(77) = since the usual view is that the effective action of massless fermion is 
non-local. One of the arguments supporting this conclusion is that the fermion operator 
In (D + 77) in free background U = I is non-local in fermionic variables in that limit, i.e. 
|| ln(D(I) + 77) n m || oc exp (—7(77) \n — m\ ) with lim^o 7(77) = 0. Indeed, for sufficiently 
small 77, the singularity of the operator in the Fourier space closest to real momenta is dis- 
tance 77 away, and hence this is also an inverse range of the operator. While it is not obvious 
a priori that the locality in fermionic and gauge variables have to be strictly related, this is 
expected to be generically true. This is expressed in the following conjecture 

Conjecture C7. Let D e S F such that f(D) = ln(D + 77) is well-defined for arbitrary 
77 > 0. If 0(77) is the inverse range of s n (U) = tr f{D{U)) n , n , then lim^ 0(77) = 0. 

8 Note that throughout the paper we could also use the definition of the massive Dirac operator which 
also scales D such that the spectrum is guaranteed to be within the fixed compact region. In case of overlap 
Dirac operator this would mean to consider (1 — rj/2p)D^ + rj instead of + 77. Such definition would 
yield all the ranges in (25) to be finite intervals. 
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Immediate issue with symmetric logarithmic LQCD (22) is that, assuming the validity 
of Conjecture C7, it is non-local by the standard definition of (strong) locality used so far 
in this discussion (see Eq. (13)). Indeed, let us analyze this in more detail. Writing the 
lattice gauge action in the standard form S (U) = (3J2 n O n (U, g) (see [1]) the requirement 
of locality is that the operator O n (U,g) satisfies Eq. (13) independently of g, at least for 
sufficiently small g. In other words, even though the optimal locality parameters A(g), a(g) 
can depend on g, we can select constants A > and a > such that (13) is satisfied 
irrespectively of g. If A(g), a(g) are continuous functions for g > 0, then this happens if 
< lim 9 _>o A(g) < oo and < lim g ^ a(g). However, in our case we have (m = wig) 

°» = U^W) tr HD + mo) - = 12^) ^ (26) 

which, according to Conjecture C7, is expected to be non-local in this strong sense because 
lim g ^ a(m (g)) = for s n- 

While strong locality is believed to be sufficient for the universality of the continuum limit, 
it is not clear that it is also necessary. Moreover, for non-ultralocal operators it is hard to 
guarantee the condition for all configurations, and one is relegated to possibly excluding a 
certain subset of them that will not satisfy it, but are not believed to be contributing in the 
continuum limit. However, the relevance in the continuum limit can in principle depend on 
the number of flavors, and possibly also on which lattice regularization is used to define it, 
which makes related arguments very hard to make rigorous. The above two points suggest 
that it is perhaps sensible to revise the requirements for locality of the theory in two ways. 
(1) Require only that the range of locality scales to zero in physical units as the continuum 
limit is approached. (2) To view the notion of locality for an operator only in conjunction 
with the theory (sequence of distribution densities) in which it is used, i.e. to assign the 
locality property to a pair (lattice operator, lattice theory) rather than to the operator alone. 

Let us now outline such definition of locality for the case of scalar (or pseudoscalar) oper- 
ators that only depend on gauge fields, since such operators are relevant for our discussion. 
According to the comments made above, the starting point in the definition involves fixing 
the lattice "theory" T considered. This, in turn, means choosing action S(p) depending on 
the set of bare lattice parameters p together with the prescription of assigning the lattice 
spacing a to the model, as well as a prescription for changing the bare parameters p = p(a,p) 
for the continuum limit a — > to be taken with some set of external physical parameters 
p being fixed. For example, in QCD with two degenerate flavors of quarks with 9 = 0, if 
we use the lattice action (22), then p = {mg,m}, the lattice spacing can be assigned e.g. 
by fixing the lowest pseudoscalar gluebal mass, and p = {frf} can be a fixed renormalized 
quark mass in some particular scheme. We can thus schematically write specification of the 
lattice theory as T(p) = { S(p), a,p(a,p) }. 

Consider an arbitrary scalar (pseudoscalar) operator O n = O n (U,p), where we allowed for 
a possible dependence of O on the same bare parameters as the theory in which its locality 
will be considered. For arbitrary configuration U let us define the function A° m (U, p) via 
(see Eq. (13)) 



dO n (U,p) 



du b (m,^) 



(27) 
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and its average in the ensemble corresponding to theory T{p) (with p being fixed and not 
explicitly denoted) at given lattice spacing a as 

A? >m (T, a) = ( A? >m (tf,p(a)) ) 5(p(a)) (28) 

where {-)s( P (a)) denotes ensemble average over the distribution specified by S(U,ip,^;p(a)) 
in the usual sense. We then consider the operator O being weakly local relative to the theory 
T{p) if the following two conditions are satisfied. 



(i) There is a lattice spacing ao > such that for all < a < ao there exist A T (a) > and 

a T (a) > such that 

A£ m (7» < A T ia)e- aT(a)ln - ml Vn,m (29) 

As before, it is implicitly assumed here that A T (a), a T (a) are chosen in an "optimal" 
way, i.e. maximal a T (a) and minimal A T (a). 

(ii) The coupling to gauge fields at arbitrary non-zero physical distance f = \n — m\a 

vanishes in the continuum limit both in absolute terms, and relative to coupling at 
f = 0, i.e. 

Urn A T ia)e- aT(a)F/a = Urn e- aT(a)p/a = Vf > (30) 

Note that condition (ii) in fact involves two requirements. In the usual situation when 
lim a ^ A T ia) = A T (0) < oo, this translates into a single requirement that the range of the 
interaction defined by a T (a) goes to zero in physical units. In other words 

lim a }. = (31) 
a^o a T ia) v ' 

In the next section we discuss the possibility that symmetric logarithmic LQCD is weakly 
local. 



4.3 Possibilities for Weak Locality 

With the acceptance of non-ultralocal operators (the conceptual leap taken decisively by 
using an overlap operator), the aspect of locality acquired a much more prominent role 
in lattice QCD than it used to. While with ultralocal operators locality was a non-issue 
decided before any investigations of dynamical behavior of the theory began, with non- 
ultralocal operators it appears that locality needs to be viewed in the dynamical context 
and, as such, is to be determined a posteriori. In this situation there exists a danger that 
one will either mistakenly put his faith in the formulation which will eventually turn out not 
to define QCD, or that one will mistakenly discard formulations with beautiful properties 
assuming improperly that they are non-local. It is thus important that one has a sensible 
"guessing guide" to make reasonable choices. It is in this context that we view the possible 
usefulness of weak locality. In other words, we propose that if lattice formulation has a 
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sensible chance to satisfy weak locality, then it also has a sensible chance to define QCD in 
the continuum limit. 

In this section we suggest that the possibility of symmetric logarithmic LQCD (22) be- 
ing weakly local depends crucially on the nature of expected divergence of c s (rj) in the 
rj — > limit. To do that, let us assume that we wish to define QCD with Nf flavors 
of quarks via symmetric logarithmic LQCD, and we will set 9 = for simplicity. We 
thus have that p = {m^ = m , mi, m 2 , • • • , m^} and, to define a theory T, we fix the 
procedure for determining the lattice spacing, as well as a scheme in which renormalized 
quark masses p = {fh[,ffi2, ■ ■ ■ ,fh r Nf } are being fixed in physical units as the continuum 
limit is taken. Conjecture C5 implies the existence of coefficients A T (a), a T (a) charac- 
terizing the operator s n (U,mo(a)) = trln(D(C/) + mo(a)) n ,n along the path to the contin- 
uum limit in this theory. We are interested in the locality properties of the gauge action 
O n (U,mo(a)) = s n (U,m Q (a))/12c s (mo(a)) (see Eq. (26)). For arbitrary p it is expected 
that A T (a) / c s (mo(a)) does not diverge for a — > since the logarithmic divergence in s n 
for m — > is removed away in O n via division by c 5 (m ). 9 At the same time the inverse 
ranges of locality for s n and O n are identical. Consequently, the validity of weak locality 
for gauge action of symmetric logarithmic QCD translates into validity of condition (31) or, 
emphasizing the implicit p dependence of the whole procedure 

l im a{ P) = (32) 

a T (m (a(p))) 

For asymptotically free theory (Nf < 16), the continuum limit a — > is realized via m — > 
due to Conjecture C6 and Eq. (21). In this regime, the localization range of s n in the 
ensemble average is expected to be crucially driven by mo irrespectively of the theory we 
are following (see the discussion introducing the Conjecture C7 as well as Conjecture C8 in 
Appendix C). In particular, a T (m ) oc m for m — > 0. Consequently, for sufficiently small 
lattice spacings, where the asymptotic 1-loop scaling formula relating a to bare coupling g 
(and hence mo) can be used, we require that 



exp(- 



\c s (m )\ 



= lim — ^— = lim m 80 (33) 

m ^0 mo m ~*0 

where fio = (H _ §iV/) /167T 2 , and the relation (21) as well as Conjecture C6 were used. 

Equation (33) implies that symmetric logarithmic LQCD for Nf asymptotically free fla- 
vors (Nf < 16) is expected to be weakly local only if 

^ > 1 (34) 

Since f3 > for Nf < 16, this means that |/t 5 (0)| > (3 which can, in turn, be viewed 
as a restriction on the number of asymptotically free flavors for which the definition via 
symmetric logarithmic LQCD is possible. In particular, 

33 

— - 24tt 2 |^(0)| < N f (35) 



9 The heuristic argument applies for smooth configurations but is not expected to be violated in the 
ensemble average. See Conjecture C8 in Appendix C. 
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We have thus arrived at a rather intriguing conclusion, namely that the suitability of 
Tr ln(D + m ) to be a gauge action with m controlling the continuum limit can in principle 
depend on the number of flavors in the theory we wish to define. In particular, there is 
possibly a minimal number of flavors for which asymptotically free SU(3) gauge theory can 
be defined in this symmetric manner. For example, pure glue QCD (Nf = 0) can only be 
formulated if 

33 

i* s (°>i > < 36 > 

It is not possible to define any asymptotically free SU(3) gauge theory via symmetric logarith- 
mic LQCD if |ft 5 (0)| < 1/487T 2 . The fact that the consistency of this maximally symmetric 
coherent LQCD might dictate the minimal number of quark flavors is an extreme example 
of how strongly the unified description of gauge and fermionic aspects of the theory can in- 
terrelate the two. Needless to say, computing k (//) for the family of overlap Dirac operators 
would shed a rather intriguing detail on this point, and this will be pursued. 
Finally, let us remark that in derivation of result (35) we have used, apart from Conjec- 
tures C4-C7, two other ingredients that are expected to manifestly hold, but are not proved. 
We give the relevant statements in the Appendix C and propose to examine their validity 
numerically for the case of overlap Dirac operator. 

4.4 Discussion 

To highlight the appeal of symmetric logarithmic LQCD, let us discuss in more detail how 
contributions of quarks and gluons to the total distribution density of gauge configurations 
appear in the completely form-symmetric manner. Consider the theory at 9 = with its 
corresponding effective action and the associated distribution density, namely 

- S tl! f} = Tr E ln { D ( U ) + m f ) = Tr In J] (D(U) + m f ) 
f=o f=o 

P {mf} oc exp(-S eff (U)) = det Y[(d(U) + m f ) (37) 

/=o 

where we used Eq. (22) and it was implicitly assumed that c s (m ) < for sufficiently 
small m . 10 Thus, in this regularization, the gauge field contribution to the total action is 
equivalent to adding an additional flavor of quarks. Indeed, one can write the total action 
of symmetric logarithmic LQCD in the form 

N f 

S{m f} = J2^ f (D(U)+m f y (38) 
f=o 

where the additional flavor ip becomes infinitely heavy (in physical units) in the contin- 
uum limit, and decouples from the light "physical" flavors in arbitrary fermionic correlation 



10 One can make very heuristic arguments for this to be the correct sign. Preliminary numerical results 
confirm this expectation at least for the overlap Dirac operator [6]. 
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functions. Indeed, sufficiently close to the continuum limit we have 



m (a) 



l- 



lim oc lim m n =00 (39) 

if Eq. (33) guaranteeing weak locality is satisfied. 

It is also instructive to emphasize the dual view, where the same lattice dynamics (at the 
effective action level) can be considered as being defined by the collection of Nf + 1 form- 
symmetric gauge action terms controlled by Nf + 1 distinct coupling constants. Indeed, the 
lattice action density can be written in the form 

S e J S = Y,-^2 F n(U;9f) (40) 
f=o z 9f 

where the (scalar) lattice operator F 2 {U;g) is defined via 

F 2 (U;g) = 2g 2 tv[\n(D{U)+r 1 {g))] nn \c s ( V )\ = ± (41) 

We note again that we assumed that c (mo) < for sufficiently small m and thus, strictly 
speaking, the above lattice formulation has the proper classical limit only for sufficiently 
small go- Among other things, form (40) nicely illustrates the point that locality should 
be viewed as a "theory-dependent" concept. Indeed, while the same operator F 2 is used to 
define the full theory, only the term driven by go (gauge action) involves weakly local operator 
along the path to the continuum limit. Indeed, for arbitrary set of fixed renormalized masses 
p = {fhi,ffi2, ■ ■ ■ ,rh r Nf }, the term driven by coupling constant gj for / > 0, will have the 
range proportional to 1/ffff in the continuum limit, and is manifestly non-local. 



5 Classically Coherent LQCD 

While our rationale for proposing coherent LQCD was based entirely on the considerations of 
QCD vacuum structure developed in Ref. [1], there is an additional and seemingly unrelated 
motivation for it. In formulating lattice regularizations we are almost exclusively guided by 
formal equations in the continuum. Indeed, we require that relevant expressions are repro- 
duced in the classical limit. Moreover, we try to arrange that symmetries of the continuum 
theory are respected by lattice dynamics to the largest extent possible. This is believed to 
ensure universality and also make the transition to the continuum limit smoother. However, 
there is another possible element in such continuum-lattice correspondence that is usually 
not taken into account. To see this, let us rewrite the continuum expression for QCD action 
(59,60) in a different manner. In particular, sorting out different Clifford components of 
D 2 = (Dp x 7 M ) 2 and D A one can easily check that 

(D 2 - D,D, x I s ) 2 i>{x) = (-\F, v {x)F, v {x) x I s + \F, v {x)F, v {x) x 7s ) ^(x) (42) 

and 

(^-(^fxl')^) = (4F P (x)F F (i) x I s + iF F (!)F F (x) x 75 + . . .) ^(x) (43) 
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One can thus formally write the continuum action as 



D^Dp x I s ) + i)(D + m)^ 

(D^D^f x I s ) + 4>(D + m)ij (44) 

Note that the above expression makes explicit the observation that we promoted in this 
article starting from very different motivation (principle of chiral ordering). In particular, 
it can be interpreted as suggesting that the fundamental object for formulation of QCD is 
massless Dirac operator D, and that gauge and fermionic aspects of the theory are explicitly 
connected to one another because of that. 

We emphasize that the above motivation for coherent LQCD is classical in nature since 
the corresponding formal equations in the continuum can be made meaningful for fields that 
are smooth almost everywhere. These equations suggest that from a classical point of view, a 
natural way to proceed with definition of coherent LQCD is to choose arbitrary lattice Dirac 
operator D e S F , and then form the coherent LQCD via the prescription given in Sees. 2 
and 3, with f(D) = D A . We will refer to this formulation as classically coherent LQCD. 
11 Note that the statement on classical limit needed for this formulation is incorporated in 
Conjecture C3g, and the locality follows from locality of D. 

We wish to make two points here. 

(i) Consider the set S F D S F of lattice Dirac operators D that satisfy all conditions specifying 
S F except the condition of lattice chiral symmetry. Then we can still follow the procedure 
above, and associate with arbitrary D e S F the scalar and pseudoscalar densities, as well as 
the gauge action via traces of f(D) = D 4 . We will refer to these gauge objects as classically 
associated with D. Thus, for example, the scalar gauge action classically associated with 
Wilson fermions is a Wilson gauge action. The use of Wilson gauge action in combination 
with overlap fermions appears highly incoherent from this point of view. Needless to say, it 
would be of interest to study the properties of the gauge action classically associated with 
overlap fermions. 

(ii) It should be emphasized here that we do not necessarily view classically coherent LQCD 
as being privileged over other choices. Rather, our view is that, at least for issues related 
to QCD vacuum structure, the selection of coherent formulation that makes the transition 
to the continuum limit smooth will depend on the nature of configurations dominating the 
QCD path integral. This is obviously an open issue. 

6 Other Operators 

Apart from large freedom in choosing the action for lattice regularization, there is an analo- 
gous freedom in choosing other lattice operators for measuring physical observables. In line 
with the point of view taken in this paper, we would like to use operators that are explicit 

11 One can also define coherent LQCD that strictly follows the form of Eq. (44). This will be discussed in 
required detail elsewhere. 



s = 



m ( 1 9 



7 5 ) {d 2 - 
75) {D* - 
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functions of lattice Dirac kernel D. This would bring the coherence (in the sense talked 
about here) to the entire process of extracting QCD predictions via lattice definition. It 
would also most likely lead to a high degree of space-time order in typical configurations of 
gauge-invariant composite fields. In fact, using the topological charge density constructed 
this way led to the basic finding that fundamental topological structure in the QCD vac- 
uum exists [11, 12]. Before we start, we should point out that all fermionic observables 
(or fermionic parts of mixed observables) are already expressed coherently since they are 
functions of (D + m) _1 , and hence D. Thus, our aim is basically to do the same for gauge 
observables. 

The generic way of achieving this is to use chiral ordering transformations of the gauge 
field defined in Ref. [1]. Indeed, to every operator 0(U) of interest, we can assign a related 
operator 

O m °{U) = 0(M D {U)) (45) 

where Ai D is the chiral ordering transformation associated with operator D. Transformation 
Ai D extracts the effective SU(3) phase acquired by chiral fermion when hopping from n + fi 
to n relative to the free case. Simple version of such transformation discussed in [1] is given 
by M D = M {3) o M {2) o M {1) with 



MW(U) = -tr s [{D ntn+ ,{l))-'D n ^{U)\ (46) 



4 

Maps M.^ and represent the unitary and group projections respectively. Clearly, the 
operator M (U) depends explicitly on D. We emphasize that for any local gauge-invariant 
operator O, the associated operator M inherits its transformation properties, locality, and 
its classical limit. Similarly, other observables of interest, such as large Wilson loops, will 
inherit their required properties. 

While the above construction achieves the goal of complete coherence for all aspects of 
the theory, it might be useful to have explicit representation for certain relevant operators 
in the same way we have for scalar density and pseudoscalar density. To do that, let us 
decompose the lattice Dirac operator D G S F in its Clifford components, namely 

16 1 

^ = ES*xr %, m = -M s rD n , m (47) 

i=i ^ 

where T = {P, % — 1, . . . , 16} is the complete orthogonal Clifford basis such that (P) 2 = P. 
We will use T = {I s , 7^, ev^), 27^75, 75}, where = ^-[7^, 7^], 75 = 7i727374, and 
7-matrices satisfy {7^,7^} = 25^ U I S . We then write the Clifford decomposition in the form 

D = §xI + V (1 X7 (1 + T F X(T p + A (1 x (^7^75) + 3> x 7s (48) 

where the gauge covariant operators {S, V M , T^-A^, 7} carry space-color indices, and have 
well-defined transformation properties under the hypercubic group. If 7-matrices are Her- 
mitian (which we assume in what follows) then §, T^, 7 are also Hermitian while V M , 
are anti-Hermitian due to 75-Hermiticity of D. While there are several possibilities for ex- 
tracting useful gauge operators from the above Clifford components, we focus here on local 
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operators obtained from D njTl . This work is already based on the fact that useful operators 
of this type can be obtained from S and 7, namely that (see Conjecture C3) 

4tr §(u(a)) oQ = C - c s a 4 tr F M;/ (0)F^(0) + £>(a 6 ) (49) 

and (see Eq. (3)) 

Atr?{u(a)) oo = -c p a 4 trF^(0)F^(0) + 0(a 6 ) (50) 

where C is a constant and U(a) is a gauge configuration obtained from arbitrary smooth 
gauge potentials A^(x) at classical lattice spacing a using prescription (2). To obtain the 
predictions for classical limits of other Clifford components, we use the same logic that is 
behind the above results. In particular, we try to identify the most general gauge covariant 
continuum functions of the gauge field with given engineering dimension and with definite 
transformation properties under hypercubic subgroup of 0(4). Such procedure proceeds by 
first finding the most general combinations of monomials in covariant derivative D p at given 
degree (engineering dimension), that do not involve any derivatives of the vector acted upon. 
The result of this straightforward manipulation obviously yields a constant for dimension 0, 
and that there are no operators of dimension 1. Moreover, linear combinations of F^ u are the 
only possibility for dimension 2, while linear combinations of [F^ u , D p ] = [F^, A p ] — d p F pv 
are the only possibility for dimension 3. This information is sufficient to deduce the leading 
terms in the classical expansion for V p , 7 pu and A p . In particular, 

4V„(l/(a)) 0)0 = -c v a 3 ([F pu ,A l/ ]-d u F pi/ ) x=Q + 0(a 5 ) (51) 

4V(f/(a)) 0)0 = tc T a 2 F pi/ (0) + 0(a 4 ) (52) 

4Ap(u(aj) ofi = -^a^e^^F^A^-d.F^)^ + 0(a 5 ) (53) 

Thus, while the operators associated with V p , A p do not appear particularly significant, 
the field-strength tensor is certainly relevant for many applications including the study of 
QCD vacuum structure. We point out that the relation equivalent to (52) has previously 
been mentioned in Ref. [13]. It needs to be emphasized that the above arguments leading 
to Eqs. (51-53) do not constitute the proof of their validity. Indeed, they rather suggest the 
conjectures analogous to Conjecture C3. Below we state this explicitly for the relevant case 
of 7 pi/ . Its validity for overlap Dirac operator will be examined in Refs. [6, 14]. 

Conjecture C9. Let A p (x) be arbitrary smooth su(3) gauge potentials on M 4 . If U(a) is 
the transcription of this field to the hypercubic lattice with classical lattice spacing a then 

tr s <vA),o(£/(a)) = ic T a 2 F pu (Q) + 0(a 4 ) (54) 

for generic D G S F . Here c T is a non-zero constant independent of A p (x) at fixed D. 

Finally, we point out that the conclusions on classical limits for Clifford components 
of D G S F described here are also expected to be valid for more general operators f(D), 
including all the cases discussed in this manuscript. 
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7 Effective LQCD II. 



Considerations on coherent LQCD offer a different viewpoint on the notion of effective LQCD 
at given fermionic response scale A F . The first form of effective LQCD, described in Ref. [1], 
is based on the eigenmode expansion of chirally ordered gauge field. 12 The basic logic of 
the construction is as follows. Starting from some original theory S(a) at lattice spacing a, 
we consider the theory S (a) obtained by transforming the ensemble of S(a) via chiral 
ordering transformation M. D (such as one based on Eq. (46)). Here D G S F is the Dirac 
operator used in S, and a(a) ~ a. The ensemble corresponding to effective theory S/Q is 
then obtained by performing this transformation with J\A D ' aAF , based on D aAp rather than 
D. Here D Xp represents an eigenmode expansion of D including eigenvalues with magnitudes 
up to Xp in lattice units. Upon taking the continuum limit, the theory defined by SjQ is 
naively expected to be described by non-local interaction with the range of locality related 
to l/Ap, but the lattice action itself is only known implicitly via its numerical ensemble. 13 
If one starts with coherent LQCD, then a natural possibility opens itself immediately for 
explicit definition of a non-local theory that could play an analogous role. Indeed, consider 
the action (6) for simplest version of coherent LQCD. At 9 = and with Nf degenerate 
quark flavors we can write the corresponding probability distribution of gauge fields as 

p K e Tr[N f MD+m)-(3D] = det [( D + m) N f e~^ D ] = JJ (A + m) Nf e~^ x (55) 

A 

If the "mother theory" is at lattice spacing a, then we define the associated effective LQCD 
at fermionic response scale A^ via 

Pa» oc n (A + m)"' e^ A (56) 

|A|<aA F 

An analogous definition is obviously possible for arbitrary coherent LQCD since the total 
action can always be eigenmode-expanded in this case. 

We finally wish to make the following two comments. 

(i) lt should be noted that it is not at all a priori obvious to what extent are effective LQCD 
I and effective LQCD II connected to one another. However, it is quite pleasing that the 
basic concept we are aiming at with effective LQCD can possibly be defined at the action 
level. 

(ii) While not entirely straightforward, it is not unreasonable to believe that the effective 
theory defined above could be directly simulated. It would be quite intriguing indeed to find 
out what kind of typical configurations would appear and how they compare to configurations 
obtained in effective LQCD I. 



12 The construction in [1] has its roots in Ref. [15, 16], and bears some technical similarities to "Laplacian 
filtering" proposed in Ref. [17] (see also [18]). 

13 We emphasize that, strictly speaking, the continuum limit at fixed hp does not necessarily have to exist. 
However, what is important for the concept of effective LQCD to be viable is that the continuum limit at 
the associated fixed momentum scale exists. These issues will be discussed in the third paper of this series. 
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8 Conclusions 



In this work we have presented a novel point of view at the process of constructing lattice 
regularizations. Considerations of universality make it possible to choose lattice dynamics 
from a very large set of possibilities (actions). This fact offers the advantage in that we 
can select the lattice theory that suits a particular problem and/or computational resources 
available. For example, large-scale lattice simulations are in majority of cases performed 
with simplest actions since it is judged that, given that they are amenable to fast simulation 
and define the correct theory in the continuum limit, it is efficient to use them even though 
they might have large cutoff effects. Taking the advantage of universality in this way, one 
typically treats the gauge and fermionic parts of the action independently of one another. 

The main message of this article comes down to the suggestion that paying attention 
to the coherence between gauge and fermionic parts of lattice dynamics might be beneficial 
in certain circumstances. More precisely, we propose that the chirally symmetric Dirac 
operator D be the unifying element in the construction of lattice actions. This conclusion 
is based on two different motivations. (1) If one accepts that the physical content of gauge 
configuration U is closely tied (or identical) to the set of effective phases affecting chiral 
fermion when hopping from n+fi to n (principle of chiral ordering [1]), then the corresponding 
transformation U — > M. D (U) drives lattice theory to the form where it is described by D. It 
is expected that configurations dominating the path integral in such theory will exhibit an 
increased degree of space-time order. (2) The formal definition for action in the continuum 
can be viewed as a matrix expression built entirely from D (see Eq. (44)) suggesting that, 
for construction of lattice actions, D can be considered a fundamental object. 

We have proposed several formulations of coherent LQCD that respect this implied co- 
herence, with the expectation that their use will make the transition to the continuum limit 
smoother at least for questions related to QCD vacuum structure. The simplest of these 
(such as theory described in Sec. 2, where S G oc (Tr D + const)) do not appear to pose ma- 
jor qualitative problems in terms of their inclusion into existing ways of simulating overlap 
fermions [7]. While this remains to be seen, it is quite clear already that there is a conceptual 
value in these formulations since they can make explicit certain aspects of QCD dynamics 
that are otherwise masked in incoherent formulations. The natural direction in this regard 
is to explore possible connections between gauge and fermionic aspects of the theory. This 
route was followed to some extent in this paper via proposing a logarithmic LQCD (Sec. 3), 
and especially symmetric logarithmic LQCD (Sec. 4). Here quarks and gluons contribute to 
the overall dynamics in a completely form-symmetric manner as expressed most clearly by 
Eq. (38). Indeed, the gluonic contribution can be viewed exactly as that of a quark that 
becomes infinitely heavy in the continuum limit. The issues related to locality properties of 
this theory were discussed in detail, and the way of resolving them was proposed. 

We have argued that, in addition to the action defining the lattice theory, one can also 
build coherently (based on D) all the operators of interest in QCD. This can be done in 
a generic way via the use of chiral ordering transformations. Moreover, the diagonal parts 
of various Clifford components associated with D offer explicit expressions for some of the 
useful composite fields (most notably in addition to scalar and pseudoscalar densities). 
Finally, a novel approach to definition of effective LQCD [1] at fermionic response scale 
was proposed. This relies on the fact that, in a coherent formulation, entire theory can 
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typically be eigenmode-expanded, and one can thus define the effective theory explicitly at 
the action level. 
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from correspondence with Mike Creutz and Tony Kennedy. This work was supported by the 
U.S. Department of Energy under the grant DE-FG05-84ER40154. 
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A Conventions in the Continuum 

Here we summarize the conventions that fix the formal equations of QCD in the continuum. 
Gauge field A p (x) e su(3) is the vector field of traceless anti-Hermitian matrices. The 
associated field-strength tensor is given by 

F^(x) = d^A v {x) - d^x) + [A lx {x),A v {x)\ (57) 

while the covariant derivative acts via 

D, <f>(x) = (<9 M + A,(x)) <f>(x) [D„ D v \ <f>(x) = F^(x) <j>(x) (58) 

The gauge part of the full action is defined by 

f 1 iO 

S G = J d i x[-—tiF^(x)F^(x) + —trF^F^x)] (59) 

where F pu (x) = \e pvpa F pcr (x), while the fermionic part for single flavor of quarks reads 

S F = J d 4 xip(x) (£>„ x lp + m) ^(x) (60) 

The above actions are invariant under local gauge transformations that take the form 

A^x) — > G{x)A p {x)G-\x) + G{x){d p G-\x)) 
tj)(x) — ► (G(i)xl s )^(x) (61) 

where G(x) G SU(3) specifies the transformation. 
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B Fermionic Determinants 



In this appendix, we emphasize a straightforward but noteworthy consequence of Conjecture 
C4. By construction, the gauge action of a given configuration in logarithmic LQCD is 
directly related to the fermionic determinant at a particular mass. Let us now apply the 
reverse logic and express the fermionic determinant for a smooth configuration in terms of 
its F 2 . Thus, consider a formal expression for the fermionic determinant in the continuum 

det (p cont (A) + to) (62) 

where A = {A^x)} are fixed smooth gauge potentials on a finite torus of size L p , and to is 
a fixed number. If we give meaning to the determinant by regularizing it on the hypercubic 
lattice using Dirac operator D e S F then we have 

det(D(A, L)+m) = e ^lHD(A,L) +m )-HD(o,L) +m) ] det (p(0, L) + m ) (63) 

where D(A,L) = D{U{ai)) represents a Dirac operator on the lattice with L sites in each 
direction, and U{ai) is a discretization of A according to prescription (2) with = L p /L. 
Note that D(0,L) denotes the free lattice Dirac operator on this lattice. If we apply the 
Conjecture C4 in the above relation we obtain 



det (D con \A) + m) = lim det(D(A, L) + 5 

= lim e~ cS{m) 4 trF^(a L n)F^(a L n)+0(al) ^ ^ + \ 

= e -c s (m)fd*xtrF^(x)F^(x) det ^ D cont^ +m ^ ^ 

Note that the above relation is still formal because the determinant of the free operator is 
not necessarily finite. However, what is physically relevant is the ratio of the determinants 
for two configurations A^ and A&\ where the free- field factor drops out at arbitrary L. We 
thus have 



det (g-%4W) + m) = s^j^^^,,^^ 



det (D cont (A^) 



m 



(65) 



where Ff£) and Ff?) are the field-strengths corresponding to A^ and A^. 

Now, instead of smooth gauge fields, consider arbitrary classical fields (i.e. with singu- 
larities allowed on the subset of space-time with measure zero), such that the corresponding 
F 2 is Riemann-integrable over the torus in question. Then the required classical limits are 
still expected to exist (see point (i) in Sec. 2) and equation (65) will still be valid. Thus, 
regardless of mass to, the relative weight of two classical configurations in QCD path integral 
is completely determined by their F 2 content if Conjecture C4 is valid. 
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C Few Relevant Statements 



For completeness (and because they are interesting in their own right) we state here a 
conjecture containing three ingredients that were implicitly used in the derivation of result 
(35). These statements relate to the existence and robust behavior of locality parameters 
A T (a), a T (a) associated with the operator s n (U,m ) = trln(D(C/) + m ). Thus, for the 
conjecture below we assume that s n (U,mo) is constructed using arbitrary D e S F such that 
\n(D + mo) is well-defined for any itlq > 0. Furthermore, we will consider the lattice action 
S(p) of symmetric logarithmic LQCD (22) for 9 = 0, in which case the set of bare lattice 
parameters consists of positive masses p = {m , m 1 , m 2 , . . . , m^ f } for arbitrary Nf > 0. 
Let us consider continuous paths p(t) (t > 0) in this parameter space that approach p = 
{0, 0, . . . , 0} monotonically in each component as t — > 0, i.e. mj(0) = and rrij(ti) < rrijfo) 
if t\ < t2, for all masses rrij. We will refer to such paths as monotonia paths for short. For 
arbitrary monotonic path we define the function A* (i) corresponding to operator s n in the 
same way as in Eqs. (27,28), namely 



K,m(t) - {K, m (U,m (t))) s{p{t)) (66) 
We propose the following conjecture to be true 

Conjecture C8. Let p(t) be a monotonic path in parameter space of symmetric logarithmic 
LQCD with arbitrary number of flavors. Then the following statements hold. 

(i) For arbitrary t > there exist positive numbers A(t), a(t) such that 

A-, m (t) < A{t)e-^ n ~ m \ Vn,m (67) 

(ii) If A{t) is optimal then the limit below exists and satisfies 

A(t) . x 

< lim , _ . v ; .. , < oo 68 
*-*> \c s (m (t))\ 

(Hi) If a(t) is optimal then the limit below exists and satisfies 

< lim < oo (69) 

t-o m (t) 



As before, by optimal pair A, a we mean maximal possible a and then minimal A given 
that choice. Note that the part (i) of the above is a consequence of Conjecture C5 (strong 
locality of s n at arbitrary fixed mo). Part (ii) asserts that the leading divergence in ensemble 
average of trln(D +m ) n ,n for m — > can be removed via division by diverging c 5 (m ). This 
conclusion is in fact stronger than what is needed for derivation of (35). The meaning of part 
(iii) is that the inverse range of locality for s n in the vicinity of m = is crucially driven 
by m not only for smooth configurations but also in the ensemble averages of symmetric 
logarithmic LQCD. 
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